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Abstract 

Background: The excessive accumulation of body fat is a major risk factor to develop a variety of metabolic 
diseases. To investigate the systematic association between the differences in gene expression profiling and 
adipose deposition, we used pig as a model, and measured the gene expression profiling of six variant adipose 
tissues in male and females from three pig breeds which display distinct fat level. 

Results: We identified various differential expressed genes among breeds, tissues and between sexes, and further 
used a clustering method to identify sets of functionally co-expression genes linked to different obesity-related 
phenotypes. Our results reveal that the subcutaneous adipose tissues mainly modulate metabolic indicators, 
nonetheless, the visceral adipose tissues as well as the intermuscular adipose tissue were mainly associated with the 
impaired inflammatory and immune response. 

Conclusions: The present study provided the evidence of gene expression profiling that the subcutaneous adipose 
tissues are mainly affected the metabolism process, whereas the visceral and intermuscular adipose tissues should 
been term as the metabolic risk factors of obesity. 
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Background 

Obesity is becoming a major global health hazard mainly 
because of the increased ingestion of high energy food 
and the decreased amounts of physical activity in popula- 
tions. A gain in the adipose mass in the body is generally 
accompanied by increased risk of metabolism-related 
diseases. However, it is not only the amount of body fat 
mass, but also the distribution of adipose tissue that is 
now recognised as an important predictor of metabolic 
abnormalities. 

Based on their anatomical location, adipose tissue can 
be divided into two major types: visceral adipose tissue 
(VAT) and subcutaneous adipose tissue (SAT). The VAT 
is a higher risk factor for obesity than SAT [1]. In par- 
ticular, basal free fatty acid flux, lipolysis rate, and 
secretory protein expression are markedly higher in vis- 
ceral adipocytes compared with subcutaneous adipocytes 
[2]. Our previous reports have provided epigenetic evi- 
dence that VAT is mainly associated with impaired 
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inflammatory and immune responses and is a metabolic 
risk factor of obesity, while SAT is mainly associated 
with metabolic processes [3,4]. In addition, ethnic group 
and gender are important factors that affect adipose 
deposition [5,6]. 

Pig (Sus scrofa) has been recognised as an attractive 
biomedical model for human because its physiology and 
anatomy are remarkably similar to humans and pigs are 
easy to breed and handle [7]. Here, we present a com- 
prehensive genome-wide comparison of gene expression 
profiling among six adipose tissues from different body 
sites in both sexes of three well-defined pig breeds 
displaying distinct adipose phenotypes. Using microarray 
technology, we identified various differentially expressed 
genes that may be related to the different adipose depos- 
ition phenotypes among the breeds, and between the 
sexes and anatomic locations. We also used a clustering 
method to identify sets of co-expressed genes that may 
be functionally linked to different obesity-related pheno- 
types. These results can be used to understand the dif- 
ferent accumulation mechanisms of the distinct adipose 
tissue types. 
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Results and discussion 

Obesity-related phenotypes 

Our previous report based on same individuals have 
demonstrated that the adipose volume were significant 
difference among adipose tissues and breeds, and be- 
tween sexes [3]. In addition, the 24 kinds of reprehensive 
metabolism indicators in serum also revealed the same 
ranking from the leaner Landrace, the wild Tibetan and 
the fatty Rongchang pigs [3] . 

Identification of differentially expressed (DE) genes and 
relatedness of adipose samples 

The number of DE genes (P < 0.05, three-way repeated- 
measures analysis of variance [ANOVA]) was lowest be- 
tween male and female (1,225), higher among different 
adipose tissues (2,750), and highest among the different 
breeds (2,949), indicating significant biological differ- 
ences between the groups in the latter two categories 
(Figure 1A). To further characterise variability in the 
gene expression profiling, we performed a principal 
components analysis (PCA) based on the DE genes in 
each of the compared groups. The six adipose tissues 
clustered into two groups, VATs and SATs, which 
reflected their distinct functional and metabolic features 
(Figure IB). Interestingly, the intermuscular adipose 
(IAD) tissue that is deposited between muscle bundles, 
clustered with the VATs in terms of gene expression 



profiling, suggesting that IAD may be a potential new 
risk factor for obesity-related diseases, in agreement with 
our findings from a previous study on methylation [4] . 

Although the three pig breeds clustered into three 
groups, the Rongchang and Tibetan breeds were closer 
to each other than the Landrace breed (Figure 1C). 
These results not only reflect the marked phenotypic dif- 
ferences in adipose deposition among the three pig 
breeds which have undergone different/opposite breed- 
ing directions, but also highlight the deep phylogenetic 
split between the European (Landrace) and Asian pigs 
(Rongchang and Tibetan) [3]. 

In addition, the male and female were also been 
roughly distinguish (Figure ID), which consisted with 
the lowest variation of gene expression profiling between 
sexes, which only have less than half of the amount of 
DE genes than among tissue types and breeds. 

DE genes between the VATs and SATs 

The PCA classified six adipose tissues into two groups: the 
VAT group which included greater omentum adipose 
(GOM), mesenteric adipose (MAD), retroperitoneal adipose 
(RAD) and IAD; and the SAT group which included upper 
layer of backfat (ULB) and inner layer of backfat (ILB). 

We performed a pairwise comparison analysis and identi- 
fied 192 and 181 DE genes that were specifically up- 
regulated in the SAT and VAT groups, respectively (Figure 2 
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Figure 1 Identification and the PCA of the DE genes. (A) DE genes across breeds, sexes, and tissues; Two-way PCA of (B) six adipose tissues; 
(C) three pig breeds; (D) sexes. L, R and T refer to the Landrace, Rongchang and Tibetan, respectively. 
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and Table 1). The DE genes which were specifically up- 
regulated in the VAT group were mainly enriched in the 
GO categories of the immune processes, including positive 
regulation of defence response (8 genes, P=5.63xl0~ 4 ), 
positive regulation of immune system process (10 genes, 
P= 1.11 x 10~ 4 ), adaptive immune response (6 genes, P = 
3.54 x 10~ 4 ), innate immune response (7 genes, P = 7.86 x 
10~ 4 ) and B cell mediated immunity (5 genes, P = 9.91 x 
10~ 4 ) (Figure 3). Various genes known to be involved in 
the immune response were highly expressed in VATs 
compared with in SATs (P < 0.05). For example, C3 (com- 
plement 3), which was reported to play a central role in 
the activation of the complement system [8,9] was up- 
regulated in the VATs. CLEC12B (C-type lectin domain 
family 12B) limits the activity of monocyte-derived im- 
mune cells after cell differentiation and possibly during 
inflammatory diseases [10] (Figure 4). The action of these 
genes which were highly expressed in VATs may contrib- 
ute to the obesity-induced chronic inflammation in adi- 
pose tissue that precedes the development of insulin 
resistance and type II diabetes. 

The DE genes that were specifically up-regulated in 
the SAT group were mainly enriched in the GO 

Table 1 Identification of DE genes which specifically 



up-regulated in SAT and VAT groups 
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Up-regulated genes 


Down-regulated genes 
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304 


291 
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154 


163 


ILB VS. GOM 


286 
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ILB VS. MAD 


260 


222 


ILB VS. RAD 


73 


60 


ILB VS. IAD 


84 
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Specific expression genes 


192 (>4) * 


181 (>4) * 



*: if a DE gene were identified from at least four comparisons, which been 
termed as the tissue- specific DE gene. 



categories of metabolic processes, including glucose meta- 
bolic process (6 genes, P= 1.16 x 10 2 ), hexose metabolic 
process (6 genes, P = 2.81 x 10" 2 ), fatty acid metabolic 
process (6 genes, P= 3.15 x 10" 2 ) and monosaccharide 
metabolic process (6 genes, P = 4.77xl0" 2 ) (Figure 3). 
This result is consistent with previous reports that, com- 
pared with the high-risk of visceral obesity caused by 
VATs, SATs are mainly associated with metabolic pro- 
cesses [4,11]. 

DE genes among the three pig breeds 

The number of DE genes (P < 0.05, three-way repeated- 
measures ANOVA) was highest between the fatty 
Rongchang and leaner Landrace pigs (1,855), lower be- 
tween the wild Tibetan and Landrace pigs (1,440), and 
lowest between the Tibetan and Rongchang pigs (868) 
(Figure 2). 

The highly expressed genes in the Landrace pigs were 
mainly involved in the GO and KEGG categories of 
skeletal system development processes and immune pro- 
cesses, including positive regulation of skeletal system 
development (7 genes, P = 2.22 x 10~ 4 ), bone develop- 
ment (8 genes, P = 4.16 x 10" 2 ), positive regulation of 
immune system process (20 genes, P= 1.19 x 10~ 3 ) and 
adaptive immune response (8 genes, P = 2.39 x 10 2 ) 
(Figure 5). The highly expressed genes in the Rongchang 
pigs were mainly involved in the processes of lipid me- 
tabolism and immune processes, such as lipid oxidation 
(8 genes, P = 1.40 x 10~ 3 ), fatty acid oxidation (8 genes, 
P = 2.53 x 10" 3 ), fatty acid metabolic process (21 genes, 
P = 4.71 x 10" 4 ), activation of immune response (6 genes, P = 
1.30 x 10" 2 ) and positive regulation of immune response 
(8 genes, P = 5.93 x 10" 3 ). This result reflects the pheno- 
typic divergence of adipose deposition between these two 
pig breeds that have undergone artificial selection for dif- 
ferent breeding targets. The Landrace breed is under 
strong artificial selection for more muscle production, 
while the Rongchang breed is selected for more adipose 
deposition. 

DE genes between the sexes 

Almost three-quarters of the DE genes between the 
sexes (879 of 1,225, 71.76%) were up-regulated in the 
males (Figure 2). These highly expressed genes were 
enriched mainly in the GO categories of the immune 
processes (Figure 6), including positive regulation of im- 
mune system process (28 genes, P = 4.49 x 10" 6 ), adaptive 
immune response (11 genes, P = 1.71 x 10~ 3 ), humoral 
immune response (7 genes, P = 1.86 x 10" 3 ), immuno- 
globulin mediated immune response (9 genes, P = 2.09 x 
10" 3 ), regulation of acute inflammatory response (5 
genes, P=1.15xl0" 2 ) and immune effector process (13 
genes, P= 1.33 x 10" 2 ). The DE genes (346 of 1,225, 
28.24%), which were up-regulated in the females 
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(Figure 2), were mainly involved in metabolic processes 
(Figure 6), including the coenzyme metabolic process 
(10 genes, P = 1.08 x 10~ 3 ), cofactor metabolic process 
(11 genes, P = 1.64 x 10" 3 ), oxidoreduction coenzyme 
metabolic process (5 genes, P=1.54xKT 2 ) and organic 
acid biosynthetic process (8 genes, P = 1.10 x 10 2 ). This 
result reflects a sexual dimorphism in body fat distribu- 
tion [12]. Females deposit relatively more fat in subcuta- 
neous/inguinal depots, which has beneficial effects on 
the control of body weight and metabolism [13]; on the 
other hand, males tend to deposit more fat in intra- 
abdominal/gonadal depots, which has been found to be 
related to risk for a number of diseases [14]. 

Interestingly, we identified some genes that potentially 
contribute to the sexual differences in obesity develop- 
ment. Compared with the females, the males exhibited 
higher expression levels of RXFP2 (P =1.19 x 10" 13 ) and 



CPB1 (P = 5.69 x 10" 6 ) (Figure 7), which have been 
related to gubernaculum testis development and sperm- 
atogenesis [15,16]. However, compared with the males, 
the females exhibited higher expression levels of PALP 
(P = 3.71 x 10" 5 ) and FABP7 (P = 8.42 x 10~ 9 ) (Figure 7), 
which have been associated mainly with the secretory 
endometrial glands and breast cancer [17,18]. 

Validation of gene expression changes by Q-PCR 

The expression changes of nine genes that are well- 
known to be related to adipose metabolism (ACADL, 
ADIPOQ, ADIPOR 7, CAV 3, MDH i, MDH 2, ME i, 
Insig 1 and SCAP) showed significant positive correla- 
tions between the Q-PCR results and their expressions 
determined by microarray (P < 0.05, Figure 8), which 
highlighted the high confidence of the results from the 
microarray approach. 
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Figure 4 Differences in mRNA expression levels of the C3 and CLEC12B gene between the VATs and SATs. 
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Figure 5 Top ten GO (Gene Ontology) terms enriched for the up- and down-regulated genes among the three pig breeds. (A) The 

Landrace versus Rongchang pig; (B) The Landrace versus Tibetan pig; (C) The Rongchang versus Tibetan pigs. The P values indicated the 
significance of the comparison, was calculated by Benjamini-corrected modified Fisher's exact test. 
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Figure 6 Top ten GO (Gene Ontology) terms enriched for the up- and down-regulated genes between the male and female pig 
breeds. The p values indicated the significance of the comparison, was calculated by Benja mini-corrected modified Fisher's exact test. 
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Figure 7 Differences in mRNA expression levels of the CPBI, RXFP2, PAEP and FABP7 gene between the between the females 
and males. 
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Figure 8 Validation of gene expression by Q-PCR. The Pearson correlation coefficient (r) and the corresponding significance value (P) were 
shown above the columns. 



Co-expressed gene sets linked to the obesity-related 
phenotypes 

Generally, quantitative characters (such as adipose de- 
position) in mammals are determined by a set of genes 
that cooperate in certain biological processes [11]. To 
further explore the relationships between the genes and 
obesity-related traits, we applied a clustering method to 
identify sets of functionally related genes linked to the 
phenotypic traits as reported previously [19]. 

We identified eight gene modules containing 947 
genes for ULB and twelve gene modules containing 
1,277 genes for ILB. In three typical VATs, we identified 
nine, nine and eleven gene modules containing 933, 
1,058 and 1,409 genes for the MAD, GOM and RAD tis- 
sues, respectively. For the IAD, we identified eleven 
modules containing 1,423 genes (Figure 9). 

Next, we explored the correlations between the mod- 
ules of co-expressed genes and the various obesity phe- 
notypes. Three gene modules in the SATs (i.e. ULB and 



ILB) were significantly correlated with adipocyte volume 
and the concentration of plasma lactate dehydrogenase 
in serum (false discovery rate adjusted permutation, 
P < 0.05) (Figure 10 and Table 2). The co-expressed genes 
within these modules were mainly involved in metabolic 
processes, such as fatty acid metabolic process (23 genes, 
P = 4.38 x 10" 5 ), fatty acid beta-oxidation (3 genes, P = 
2.65 x 10" 3 ), coenzyme metabolic process (6 genes, 
P = 4.77 x 10' 5 ) and carbohydrate metabolic process 
(12 genes, P = 2.94 x 10~ 3 ) (Table 3), in agreement with an 
earlier report that SATs contribute mainly to metabolic 
processes [20,21]. 

The co-expressed genes in two of the VATs (i.e. GOM 
and MAD) mainly affected the concentrations of inflam- 
matory and immune adipokines in serum, i.e. plasmino- 
gen activator inhibitor- 1 (PAI-1), interleukin 6 (IL-6) 
and tumor necrosis factor-a (TNF-a) (Figure 10 and 
Table 2). The co-expressed genes within these modules 
were mainly involved in immune system process (21 
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correlations between genes residing in all the gene sets were plotted. Bright red indicates a strong positive co-expression, whereas bright green 
indicates strong negative co-expression. 



genes, P = 6.66 x 10" ), immune response (8 genes, P = 
6.84 x 10" 3 ), cellular defence response (15 genes, P =1.25 x 
10" 10 ), cellular defence response (6 genes, 
and B cell mediated immunity (9 genes, P=1.05xl0~ 7 ). 
This result suggested that the GOM and MAD tissues con- 
tributed mainly to the obesity-induced chronic inflamma- 
tion in adipose tissue that precedes the development of 
insulin resistance and type II diabetes [22,23]. The RAD tis- 
sue, on the other hand, was found to affect mainly the body 
density and the concentrations of three metabolic indica- 
tors in serum, namely, high density lipoprotein, insulin, and 
apolipoprotein A-l (Apo-Al), which are mainly involved in 
metabolic processes such as coenzyme metabolic process 
(16 genes, P =3.02 x 10" 9 ), primary metabolic process (241 
genes, P = 9.31 x 10" 9 ), and tricarboxylic acid cycle (7 genes, 
P=5.49x 10" 6 ) (Table 3). Studies in human showed that 
GOM was strongly correlated to insulin, while no such 
correlations were found between lipid uptake and SAT or 
RAD [24]. 

Similar to the findings for the high-risk VATs, the co- 
expressed genes in IAD were also found to be significantly 



correlated with immune- related TNF-a (Figure 10 and 
Table 2), further suggesting that the IAD is an independ- 
ent risk factor for metabolic diseases [25]. 

Conclusions 

The present study describes a genome-wide analysis of 
gene expression profiling among various adipose depots 
between the sexes of three well-defined pig breeds 
displaying distinct adipose phenotypes. We present evi- 
dence that the SATs mainly affect metabolic processes, 
while the VATs and the IAD are associated mainly with 
immune and inflammation responses that are seen as 
the metabolic risk factors of obesity. This study also sup- 
ports pig as an ideal model for studying human obesity. 

Materials and methods 

Tissue collection 

Three females and three males at 210-days-old for each 
of the leaner Landrace pigs (a Western breed), the wild 
Tibetan pigs (a feral, indigenous Chinese pig that has not 
undergone artificial selection) and the fatty Rongchang 
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Modules 

Figure 10 Correlations between co-expressed gene modules in six adipose tissues and obesity phenotypic traits. -Log 10 P values for 
Spearman rank correlation coefficients between the expressed values of the gene modules and the different obesity phenotypic traits are shown. 
The gray line represents the Bonferroni corrected P value greater than 0.05. 



pigs (a Chinese breed) were used in this study as previously 
described [3]. Animals were humanely sacrificed, in com- 
pliance with experimental animals established by the Min- 
istry of Agriculture of China. Three VATs (i.e. GOM, MAD 
and RAD), two SATs (i.e. ULB and ILB) as well as IAD 
were rapidly separated from each carcass, immediately fro- 
zen in liquid nitrogen, and stored at -80°C until RNA ex- 
traction. For more information, please refer to Li et al [3]. 

Measurements of the obesity-related phenotypes 

Measurements of the pig body density, the concentra- 
tions of 24 serum- circulating indicators of metabolism 



and the adipocyte volume are from our previous report 
based on same individuals. For more information, please 
refer to Li et al [3] . 

RNA isolation and microarray analysis 

The total RNA (10 \ig) of the 108 tissue samples was 
extracted with TRIzol (Invitrogen) and further purified 
using an Rneasy column (Qiagen). Hybridization, wash- 
ing, and scanning were done according to standard 
Agilent protocols. The data analysis was performed with 
MultiExperiment Viewer (MeV) [24]. Array data have 
been uploaded to NCBIs Gene Expression Omnibus 



Table 2 Modules correlated with the obesity-related traits in SATs and VATs 



Tissues 


Module ID 


Correlated trait 


No. of genes in module 


Correlation coefficient 


P-value 


SAT ULB 


4 


Adipocyte volume 


110 


-0.85 


8.12 X10" 6 




6 


LDH 


72 


-0.84 


1 .23 x 1 0" 5 


ILB 


2 


Adipocyte volume 


157 


0.87 


2.05 x 10" 6 


VAT GOM 


2 


PAI-1 


162 


-0.81 


4.97 x 10" 5 




3 


IL-6 


92 


-0.80 


6.28 X10" 5 


MAD 


2 


TNF-a 


152 


0.82 


2.57 x 10" 5 


RAD 


8 


Body density 


71 


0.87 


2.97 x 10" 6 




11 


Body density 


51 


-0.81 


5.44 x 10" 5 




1 


HDL 


590 


0.85 


7.42 x 10" 6 




11 


Apo-A1 


51 


-0.81 


3.91 X10" 5 




3 


Insulin 


118 


0.77 


6.89 X10" 6 


IAD 


10 


Body density 


55 


0.86 


5.90 X10" 6 




11 


TNF-a 


52 


-0.83 


2.15 x 10~ 5 
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Table 3 Top ten GO (Gene Ontology) and pathway categories enriched for co-expressed gene modules which 
correlated with obesity phenotypic traits 



Correlated trait Tissue 

(Gene sets order No.) 


Functional 
category 


Term 


P values 


Involved 
gene No. 


Adipocyte volume ULB (4) 


GO-BP 


metabolic process 


1.19 x 10 7 


119 


ILB (2) 


GO-BP 


coenzyme metabolic process 


8.84 x 1CJ 6 


8 




GO-BP 


lipid metabolic process 


4.38 x 1CJ 5 


23 




GO-BP 


fatty acid metabolic process 


8.95 x 1CJ 5 


10 




GO-BP 


primary metabolic process 


1.07 x 1CJ 4 


104 




GO-BP 


translation 


5.37 x10~ 4 


12 




GO-BP 


protein folding 


1 .50 x 1 0 3 


8 




GO-BP 


fatty acid beta-oxidation 


2.65 x 10~ 3 


3 




GO-BP 


pentose-phosphate shunt 


3.51 x10~ 3 


2 




GO-BP 


generation of precursor metabolites and energy 


4.15 x 10~ 3 


9 


Body density RAD (8,11) 


GO-BP 


translation 


1.44 x 10~ 3 


8 


IAD (10) 


GO-BP 


protein targeting 


1.90X10" 3 


5 




GO-BP 


lysosomal transport 


5.44 x 10~ 3 


2 




GO-MF 


structural constituent of ribosome 


6.10 x 10~ 3 


5 




GO-MF 


structural molecule activity 


6.62 x 10~ 3 


15 




GO-BP 


nitrogen utilization 


1.20 x 10~ 2 


1 




KEGG 


Serine glycine biosynthesis 


2.97 x 10~ 2 


1 




GO-BP 


regulation of amino acid metabolic process 


3.55 x 10~ 2 


1 




GO-BP 


pyrimidine base metabolic process 


4.40 x 10~ 2 






GO-BP 


pentose-phosphate shunt 


4.71 x 10~ 2 


1 


Apo-AI RAD (11) 

VIVIcLdUUIIV. IIIUILdlLM 5) 


GO-BP 


nitrogen utilization 


3.51 x 10~ 3 


1 


GO-BP 


regulation of amino acid metabolic process 


1 .05 x 1 0 2 


1 




GO-BP 


lysosomal transport 


3.12 x 10~ 2 


1 




GO-BP 


nitrogen compound metabolic process 


3.12 x 10~ 2 


1 




GO-BP 


neuromuscular synaptic transmission 


4.64 x 10~ 2 


1 




GO-MF 


hormone activity 


1.55 x 10 ^ 






KEGG 


Vasopressin synthesis 


2.09 x 10~ 2 


1 




KEGG 


De novo pyrmidineribonucleotidesbiosythesis 


3.12 x 10 ^ 


1 




KEGG 


Cortocotropin releasing factor receptor signaling pathway 


5.14 x 10~ 2 


1 




KEGG 


De novo purine biosynthesis 


5.64 x 10~ 2 


1 


HDL (Metabolic RAD (1) 

IllUICdlUi 3) 


GO-BP 


generation of precursor metabolites and energy 


1.20X10" 15 


37 


GO-BP 


metabolic process 


1.49x10~ 14 


273 




GO-BP 


respiratory electron transport chain 


1.35x10~ 11 


30 




GO-BP 


coenzyme metabolic process 


3.02 x 10~ 9 


16 




GO-BP 


primary metabolic process 


9.31 x 10~ 9 


241 




GO-BP 


signal transduction 


6.95 x 10~ 7 


48 




GO-BP 


cell communication 


8.92 x 10~ 7 


52 




GO-BP 


cell surface receptor linked signal transduction 


1.49 x 10~ 6 


15 




GO-BP 


system development 


3.58 x 10~ 6 


9 




GO-BP 


tricarboxylic acid cycle 


5.49 x 10~ 6 


7 


LDH (Metabolic ULB (3) 


GO-BP 


fatty acid metabolic process 


1.37 x 10~ 5 


9 


indicators) 


GO-BP 


metabolic process 


1.86 x 10~ 5 


75 




GO-BP 


fatty acid beta-oxidation 


3.05 x 10~ 5 


4 




GO-BP 


coenzyme metabolic process 


4.77 x 10~ 5 


6 




GO-BP 


cellular amino acid biosynthetic process 


9.01 x 10~ 5 


5 




GO-BP 


cellular amino acid metabolic process 


1 .05 x 1 0 -4 


9 




GO-BP 


protein folding 


4.66 x 10~ 4 


7 
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Table 3 Top ten GO (Gene Ontology) and pathway categories enriched for co-expressed gene modules which 
correlated with obesity phenotypic traits (Continued) 





GO-BP 


pentose-phosphate shunt 


1.45 x 10" 3 


2 


GO-BP 


primary metabolic process 


1.99 x 1CJ 3 


66 


GO-BP 


carbohydrate metabolic process 


2.94 x 1CT 3 


12 


Insulin (Metabolic RAD (3) 
indicators 


KEGG 


VEGF signaling pathway 


3.72 x 1CJ 4 


4 


GO-BP 


transcription, RNA-dependent 


4.06 x 1(T 4 


2 


GO-BP 


cell communication 


6.62 x 1CT 4 


31 


GO-BP 


mesoderm development 


1.04 x 1CJ 3 


11 


GO-BP 


cellular process 


1.10 x 10 3 


40 


GO-BP 


protein amino acid ADP-ribosylation 


1.11 x 10~ 3 


2 


GO-BP 


signal transduction 


1.37X10" 3 


29 


GO-BP 


esicle-mediated transport 


1.69X10" 3 


11 


GO-BP 


system development 


1 .93 x 1 0" 3 


14 


GO-BP 


cell surface receptor linked signal transduction 


4.81 x 10~ 3 


16 


IL-6 (Inflammatory and GOM (3) 


GO-BP 


immune system process 


6.66 x 10~ 4 


21 


immune adipokines) )AD ^ ^ 


GO-BP 


response to stimulus 


7.72 x 10~ 4 


16 


GO-BP 


response to interferon-gamma 


2.08 x 10~ 3 


4 


GO-BP 


cellular defense response 


4.67 x 10~ 3 


6 


rr\ dd 


immune response 


O.o4 X I U 


8 


GO-BP 


DNA recombination 


1 .00 X I 0 


3 


oU-br 


nitrogen utilization 


1 .1 6 X 1 0 2 


1 


CC\ DD 


gluconeogenesis 


1 .35 X 1 0 2 


2 


GO-BP 


antigen processing and presentation of peptide or 
polysaccharide antigen via MHC class II 


i Qi v 1 rr2 
1 .0 I X I U 


2 


KEGG 


T cell activation 


1 .84 x 1 0~ 2 


4 


TNF-a (Inflammatory and MAD (2) 
immune adipokines) 


GO-BP 
GO-BP 


antigen processing and presentation of peptide or 
polysaccharide antigen via MHC class II 

antigen processing and presentation 


4.26 x 10~ 5 
6.56 x 10~ 4 


4 
4 


KEGG 


T cell activation 


1.37 x 10~ 3 


4 


GO-BP 


cellular defense response 


3.14 x 10~ 3 


6 


GO-BP 


immune system process 


8.14 x 10~ 3 


17 


GO-BP 


tRNAaminoacylation for protein translation 


2.18 x 10~ 2 


2 


KEGG 


5-Hydroxytryptamine biosynthesis 


3.17 x 10~ 2 


1 


KEGG 


Mannose metabolism 


3.69 x 10~ 2 


1 


GO-BF 


establishment or maintenance of chromatin architecture 


3.85 x10~ 2 


4 


PAI-1 (Inflammatory and G0M(2) 
immune adipokines) 


GO-BP 


immune system process 


6.00x10^ 5 


41 


GO-BP 


response to stimulus 


6.31 x 10" 15 


30 


GO-BP 


cellular defense response 


1.25 x 10 10 


15 


GO-BP 


immune response 


1.35 x10~ 10 


19 


GO-BP 


B cell mediated immunity 


1 .05 x 1 0 -7 


9 


GO-BP 


natural killer cell activation 


4.17 x 10~ 6 


7 


GO-BP 


antigen processing and presentation 


6.62 x 10~ 6 


6 


KEGG 


T cell activation 


2.04 x 10~ 5 


6 


GO-BP 


cell surface receptor linked signal transduction 


1.52 x 10^ 


23 


KEGG 


Inflammation mediated by chemokine and cytokine signaling 
pathway 


3.11 x 10~ 4 


8 


In all tests, the unified set of co-expressed genes for different adipose depots correlated to an obesity phenotypic trait was compared with the whole known 


genes which were appointed as the background. P values indicating significance of the overlap between various genes sets, was calculated using Benjamini 




corrected modified Fisher's exact test. BP biological process, MF molecular function. 
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(GEO) [accession number GSE30343]. For more infor- 
mation, please refer to Li et al [3] . 

Identification of differentially expressed genes 

We used three-way repeated-measures ANOVA (n = 3 
per breed per sex per tissue) to identified differentially 
expressed (DE) mRNA among six adipose tissues be- 
tween sexes of three pig breeds (P < 0.05, corrected with 
adjusted Bonferroni method, FDR <0.05, 1,000 permuta- 
tions), which were further analyzed using DAVID to 
examine whether these DE genes were enriched for spe- 
cific functional catalogs (i.e. KEGG pathways and Gene 
Ontology-Biological Processes and -Molecular Func- 
tions) (P < 0.05, corrected with adjusted Bonferroni 
method) [26]. 



in a set that was assigned the average expression value 
of the two genes that constitute this set. After addition 
of this newly created set to the dataset, the two individ- 
ual genes were removed from the data and the strongest 
correlation in the dataset was again selected. This 
resulted in either the expansion of a set already created 
or in the creation of a new set. We kept repeating this as 
an iterative process until the most significantly corre- 
lated pair was coefficient of correlation r < 0.80. Only the 
sets containing 50 or more genes are kept for further 
analysis. The gene sets - reflecting the average expres- 
sion value of the genes constituting that set - were cor- 
related with obesity phenotypic traits using a Spearman 
rank correlation coefficient. 



Identification of co-expressed gene modules 


Quantitative PCR 








To identify sets of functionally related genes in six adi- 


Quantitative PCR 


(Q-PCR) was 


used to 


confirm the 


pose tissues, linked to the phenotypic traits, we applied 


expression pattern observed in the microarray results. Q- 


a clustering 


method as described previously [4,19] with 


PCR was performed on a CF96 Real-Time PCR Detection 


some modifications. Spearman rank correlation coeffi- 


System 


(Bio-Rad) 


using SYBR® 


Green Real-time PCR 


cients were 


determined between all possible gene-pairs 


Master 


Mix (TaKaRa, China). 


The PCR 


primer se- 


across 18 samples for each adipose tissue. The strongest 


quences 


are shown 


in Table 4. Porcine ACTB, TBP and 


correlated g 


;ene-pair was selected, and grouped together 


TOP2B were simultaneously used 


as endogenous control 


Table 4 Primers used for the Q-PCR assay 












Gene symbol 


Primer sequence (5'— > 3') 




Amplicon length (bp) 




GenBank No. 


ACTB* 


F: TCTGGCACCACACCTOT 




114 






DQ178122 


R: TGATCTGGGTCATC1TCTCAC 


TBP* 


F: GATGGACGTOGGmAGG 




124 






DQ178129 


R: AGCAGCACAGTACGAGCAA 


TOP2B* 


F: AACTGGATGATGCTAATGATGCT 




137 






AF222921 


R: TGGAAAAACTCCGTATCTGTCTC 


ACADL 


F: GTAAGAACAAATGCCAAGA 




103 






NM_2 13897 


R: CAGCCACTACAATCACAAC 


ADIPOQ 


F: GGGTCACTGTCCCTAAC 




209 






NM_214370 


R: GTCCTGGTACTGGTCGT 


ADIPOR1 


F: TGGCTGAAGGACAACGAC 




228 






NM_001007193 


R: CAAGAAGAACATCCCAAACAC 


CAV3 


F: GCATCAGCCATATCTACTCACT 




107 






NM_001037149 


R: CAOTCmCCGCAGCAT 


MDH1 


F: TAAGGTOTCGTGGTGGG 




124 






U44846 


R: TGCmAGCTCGGTOTG 


MDH2 


F: CGAGGTGGTCAAGGCTAAG 




172 






M 16427 


R: CAATGGCGTGGAGAAATAC 


ME 1 


F: GTOCCOTGGTGTOT 




212 






X93016 


R: GGATAAATGGTGGCTGTC 


Insig 7 


F: TCAACCACGCCAGTGCTAA 




251 






AY336601.2 


R: CACGCCTCCCGAGAAGAAA 


SCAP 


F: CGCCTGAAACAGAAATCG 




213 






XM_00 1924223 


R: TCTCCTGAGCCTCCAACA 



*: ACTB (p actin), TBcP (TATA box binding protein) and TOP2B (topoisomerase II (3) are the endogenous control genes. 
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genes. The 2~ AACt method was used to determine the 
relative mRNA abundance for the surveyed samples. 

Abbreviations 

SAT: Subcutaneous adipose tissue; VAT: Visceral adipose tissue; ULB: Upper 
layer of backfat; ILB: Inner layer of backfat; GOM: Greater omentum adipose; 
MAD: Mesenteric adipose; RAD: Retroperitoneal adipose; IAD: Intermuscular 
adipose. 
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